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Critical parameters for quantum wells 2 

Abstract. We apply the Riccati-Pade method to one-dimensional quantum wells. 
The approach is based on a rational approximation to the logarithmic derivative of the 
wavefunction. We modify the Riccati-Pade method and calculate accurate critical 
parameters for the Gaussian well and the Yukawa potential. The application of 
perturbation theory to a rational potentials yields interesting and baffling unphysical 
results. 



1. Introduction 

The accurate calculation of the number of bound states supported by a finite quantum- 
mechanical potential well is of great physical and mathematical importance and for 
this reason there has been considerable interest in the derivation of upper and lower 
bounds p14T3] . Most of those bounds are given in terms of the potential-energy function. 
In a recent paper Liverts and Barnea [1^ proceeded in a different way and proposed 
the calculation of the critical parameters for negative central-field quantum wells. To 
this end they applied two exact methods and the WKB approach, the latter for the 
estimation of the large-quantum number behaviour of the critical parameters. In this 
context a critical parameter is the value of a potential parameter for which an energy 
eigenvalue is exactly zero (what the authors call a transition state). Obviously, if we 
know a critical parameter for a given potential we also know the number of bound states 
supported by it. 

Local methods, like, for example, a power-series approach, are expected to be 
unsuitable for the calculation of critical parameters. Even the Riccati-Pade method 
(RPM [l5l|T6]), based on Fade approximants was shown to be unsuccessful for the 
calculation of the eigenvalues of the Yukawa potential close to the zero-energy threshold 
(transition state) [15]. The purpose of this paper is to investigate in more detail whether 
those earlier results already prove that the RPM is actually useless for the calculation 
of critical parameters. 

In section |2] we outline the main ideas of the RPM. In section |3] we briefly discuss 
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the solutions of the Schrodinger equation with even-parity potential wells. In section H] 
we apply the approach to some simple one-dimensional models. In section [5] we propose 
a modified version of the RPM that is suitable for the calculation of critical parameters 
and apply it to the Gaussian well. In section [6] we apply the modified RPM to the s- 
states of central-field models and choose the Yukawa potential as a suitable illustrative 
example. In section [7] we discuss the application of perturbation theory to a model 
with a rational potential that is exactly solvable at threshold. Finally, in section [8] we 
summarize the main results and draw conclusions. 

2. The Riccati-Pade method 

We consider the eigenvalue equation 

ijj" {x) + Q{E,x)tlj{x) = 0, -oo < X < oo (1) 
were E is the eigenvalue. We assume that iIj{x) can be expanded about the origin as 

oo 

^(x) = x'J2 9^'^^ a,(3>0 (2) 
i=o 

It is clear that 

can be expanded about the origin as 

oo 

/(x) = ^ z = (4) 

We approximate /(x) by means of a rational function of the form x'^~^[M/N]{z) 
where 

[M/N]{z) = ^'=\' . = T{M + N + l,z) + 0(z^+^+2) 



and 



nn,z) = j:fjz^ (6) 

j=0 

We choose M > N and define d = M — N. It is not possible to satisfy the condition 
(E]) unless Hf) = \fi^j^d-i\fj=i = 0, D = N + 1. The coefficients fj, and thereby the 
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Hankel determinant H'^, depend on the eigenvalue E. From the roots of Hf){E) = we 
obtain sequences 

^[D,d]^ D = 2,3,... that converge toward the eigenvalues of Eq. ([T]) . 
The ordinary Fade approximation to f{x) is 

[M/N]{z) = ^ff^ = T{M + N,z) + 0(^*^+^+1) (7) 

If Zq is a zero of the denominator then 



^0 ~ u I • • • ~ I — iV-T \°) 

On onZq OnZq 

Suppose that E* is a zero of b]\f{E) = and that 6iv-i(-E') does not exhibit a zero 
in the interval {E* — e, E* + e). Therefore \zo\ — oo as E E*. The coefficient 
6 AT is proportional to the Hankel determinant H'^\{E) so that the Hankel condition 
Hf){E) = is equivalent to moving a pole of a rational approximation towards infinity 
[T7] . It is also equivalent to moving a zero of the approximate iplx) towards infinity. 
For this reason the Hankel condition is equivalent to selecting bound states that vanish 
at infinity. 

3. Parity-invariant finite wells 

The Riccati-Fade method is known to produce accurate eigenvalues for infinite wells or 
sufficiently deep finite ones. The purpose of this section is to investigate to which extent 
it is possible to apply the RFM to shallow wells. To this end we consider the eigenvalue 
equation ([T]) with 

Q{x) = 2[E- Vix)] (9) 

where the potential-energy function V{x) exhibits a minimum at V{0) < and 
V{x — 7- ±oo) = 0. In order to simplify the discussion we assume that V{—x) = V{x) so 
that 

oo 

Qix) = Y.Q,x'^ (10) 

j=0 

The shape of a parity-invariant potential is commonly determined by a smaller number 
of parameters. In addition to it, the results for the odd states also apply to the solutions 



Critical parameters for quantum wells 5 
of the Schrodinger equation with a central- field potential having zero angular momentum 
quantum number / (s states). 

The eigenfunctions of the Hamiltonian operator with a parity-invariant potential 
are even or odd; therefore, s = for the former, s = 1 for the latter and ^ = 2 in both 
cases. We thus have 

oo 

f{x)^xY.fjx'^ (11) 

j=0 

where the first coefficients are 



f _ Qo 

JO — 



l + 2s 



n ^ Ql , Ql 

(3 + 2s)(H-2s)2 3 + 2s 

J2 



(5 + 2s) (3 + 2s) (1 + 2s)3 (5 + 2s) (1 + 2s) (3 + 2s) (5 + 2s) 

(12) 

Besides, the function f{x) is a solution to the Riccati equation 

fix) + -fix) - fixf - Q{x) = (13) 

X 

For convenience we define vq — — V^(0) > and v{x) — —V{x)/vo. For all values 
of the well depth vq > there is always a bound ground state with energy Eq. The 
number of bound states with energies Eq < Ei < . . . < < depends on vq. As Vq 
decreases the highest bound-state energy, say E^ approaches the threshold E — from 
below and we have a critical well parameter vo,n when E^ — 0. 

We assume that V{x) vanishes faster than x~^; that is to say 

lim x^V{x) = (14) 
Therefore, an eigenfunction for the arbitrary energy E < behaves asymptotically as 

V^(x) - A{E, i;o)e-"^ + B{E, ^;o)e"^, \x\ oo (15) 

where a — \/—1E. The bound states are given by the condition B{Ej,vo) — that 
leads to square-integrable eigenfunctions. 
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When E = the two asymptotic solutions in Eq. f llSI) are hnearly dependent. In 
this case the general solution to ip'^x) = behaves as 

ipix) ^ A{vo) + B{vo)x, E = 0, \x\ ^ oo (16) 

The solution at threshold is not square integrable but we can think of it as the limit of 
a square integrable one limE^o A{E , vo)e~°'^ = A{vo) {\x\ — )■ oo). Therefore, the critical 
parameters are roots of B{vo^n) = and the boundary condition at threshold is 

lim ip'{x) = 0, E = (17) 



In their study of central-field quantum wells Liverts and Barnea [Mj assumed the 
boundary condition limj._j.oo IV'P = for the transition states with E = 0. This boundary 
condition may be valid for nonzero angular-momentum quantum number (/ > 0) but it 
does not apply to / = where the asymptotic behaviour is similar to f|T6l) . In fact their 
equation (24) for Z = is similar to our equation (|T71) . 

4. Examples 

In what follows we discuss some simple model potentials to illustrate the application of 
the RPM. 

4.I. Modified Poschl- Teller potential 

As a first example we consider the modified Poschl- Teller potential 

V{x) = ^ vo > (18) 

cosh [X) 

There are two reason for this choice: first, we can solve the Schrodinger equation and 
obtain a simple expression for the eigenvalues: [T8] 

En= ^, n = 0,l,... 

A =l±^/^ (19) 

Second, the exact bound-state solutions are hypergeometric functions of y = cosh^ x so 
that the RPM based on an x-power series can only yield approximate results. Therefore, 
this model is a suitable benchmark for testing the performance of the approach. 
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As discussed above the critical values of the potential parameter vq are determined 
by the condition i?„(t>o,n) = 0. It follows from equation ( fT9l) that fo,n = n{n + l)/2, 



Hankel determinant Hf) as vq approaches fo,o = and fo,2 = 3. We calculated E^^''^\ 
D = 2, 3, . . . for Vq = fo,n + lO^'^, for k = 1,2,.... The results show that there are 
convergent sequences of roots for D = Dk,Dk + 1, . . . and that the starting point of 
each sequence Dk increases as k increases. There seems to be convergent sequences no 
matter how large the value of k. In other words, the RPM appears to be successful no 
matter how close vq is to the critical value vo,n- 

Since the roots of the Hankel determinants H'^{E,vo) = give rise to sequences 
that clearly converge towards the eigenvalues -E„(fo) for vq quite close to fo,„, then one 
would expect to find sequences of roots of Hf){E = 0,vo) = that converge towards 
the critical parameters fo,n- This is exactly the case for this model and one obtains 
the critical parameters with any desired degree of accuracy with Hankel determinants 
of relatively small dimension. There are, however, two surprising facts. The first one is 
that the RPM yields all the critical parameters fo,„, n = 1,2,... when choosing either 
the even (s = 0) or odd {s = 1) functions. One would expect the even or odd values of n 
to appear separately with even or odd functions, respectively. The second surprising fact 
is that the RPM with s = yields the critical parameters with odd n more accurately 
than those with even n. The opposite situation takes place when choosing s = 1. 

We can understand the occurrence of twice as much critical parameters as expected 
by obtaining the corresponding wave functions in the usual way [18]. If 'ipn,s{x) denotes 
the solution of parity s for vq = wo,n and E = 0, then the first of them are given by 



n = 0, 1, . . .. We first investigate if there are converging sequences of roots El-^'^l of the 




2 



+ 1 

2 (4e2^ - e^^ - 1) 
g4x ^ 2e2^ + 1 
[e^^ (2x - 3) - 8a;e2^ + 2a; + 3] 
4 (e^^ + 2e2^ + 1) 



(20) 
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We appreciate that ipi^i{x) and 4'2fi{x) are convergent while ipifii^) and "^2,1(3^) are 
divergent. In general, ilJn,s is convergent or divergent provided that n + s is even or odd, 
respectively: 

lim ilJn,sX^) = A, n + s = 2k 
lim x~'^tpn,s{x) =B,n + s = 2k — 1 

\x\^oo 

k =1,2,... (21) 

We conclude that the RPM approaches both the convergent and divergent solutions for 
this problem when E = 0. This is the reason why the whole set of critical parameters 
vo,n appears for both the even and odd solutions: half of them are convergent and the 
other half divergent. It is clear, as already argued above, that the RPM does not yield 
the exact result because f{x) is not a rational function of x for any of the functions ( l20l) . 

It is not clear to us why the RPM with s = (s = 1) yields the critical parameters 
with odd (even) n more accurately. We will discuss this point with somewhat more 
detail later on by means of another solvable model. 

4-2. Gaussian well 
The Gaussian well 

Vix) = -v,e-^" (22) 

is another suitable choice because the potential is extremely simple but the Schrodinger 
equation is not exactly solvable. In this case the behaviour of the sequences of roots 
of Hf){E,VQ) = appears to be similar except that the starting point of a given 
sequence increases more pronouncedly as vq approaches fo,n and we could not find 
converging sequences of roots of Hfj[E = 0,fo) = 0. It is not clear to us which is 
the feature of this well that makes such a difference. The present form of the RPM is 
unsuitable for obtaining the critical parameters for this problem and in Sec. [5] we will 
discuss an improved version of the approach. 
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4-3. Rational potential 

The third example in this section is the potential well 

V(x) = (23) 

that satisfies the condition f lT4|) . Joseph [19] studied the family of central-field potentials 
V{r) = — Ar"~^(rQ + r^)~" in his discussion of local degeneracy. Clearly, the potential 
(|23|) is the one-dimensional version of the case a = 2. Besides, present results for the 
odd states should agree with those obtained by Joseph for a = 2 and / = 0. 

The roots of the Hankel determinants H%{E = 0, vq) yield exact critical parameters 
Vo^n = n{n + 2)/2, n = 1,2, . . .. These results correspond to exact rational solutions 
to the Riccati equation f ll3p . In order to understand their meaning we construct the 
corresponding wavefunctions as 



V'n,s = exp 
The first even and odd ones are 



f{x) dx 



(24) 



^i,o(a;) 
^^2,0(2;) 
^^3,0(2;) 



1-3x2 



X2 + 1 

(x^ + 2x - 1) (x^ - 2x - 1) 

(1 + x2)^/2 



5x^ - lOx^ + 1 

V'4,o(a;) = — — — — (25) 

1 + X"' 



and 



^1,1 (a;) 



X 





x^ 


X (x^ 


-3) 


1 + 


x^ 


X (x^ 


-1) 



:i+x2)'/' 



xix^ — lOx^ 5) , , 

^4,1 (x) = ^ ^ ^ (26) 

1 + x^ 



Critical parameters for quantum wells 10 
respectively. As in the case of the modified Poschl- Teller potential the solutions ipn,s 
are convergent or divergent provided that n + s is even or odd, respectively; more 
precisely, they satisfy equations (pTj) . According to the discussion of Sec. ^i'2k,o{x) are 
the even solutions to the Schrodinger equation ioi E = and vq = 2k{k + 1) (satisfy the 
condition B{vo) = 0). On the other hand, 4'2k-i,i{x) are the odd solutions for E = and 
Vq = (4fc^ — l)/2, k = 1,2, . . .. The latter agree with Joseph's ones when A = 2vo |19j . 
The remaining solutions ip2k,i{x) and ip2k-i,o{x) are the unphysical divergent solutions 
to the Schrodinger equation. We see that the RPM yields the exact convergent and 
divergent solutions to the Schrodinger equation with the potential ( l23l) when E = 0. 
It is worth noting that the RPM does not distinguish between physical an unphysical 
results unless one manages to obtain the wavefunction from its logarithmic derivative 
as we did it in this example. In other cases, like the potential ( IT8l) . it may be easier 
to resort to another approach to obtain the wavefunction from the parameters given by 
the RPM. 

In the appendix we solve the Schrodinger equation for this potential and derive the 
exact convergent and divergent solutions for E = 0. 

Although the RPM yields the exact critical parameters it is not suitable for the 
calculation of the energies close to threshold. The sequence of roots of H'^{E,vo) = 
converge rather too slowly when vq is close (slightly greater than) a critical parameter. 
When E < the function f{x) is not an exact rational function and the RPM yields 
approximately those eigenvalues that are not too close to threshold. We calculated the 
ground state for Vq = 3/2 {Ei = 0) and the first two bound states for f o = 4 {E2 = 0). In 
the three cases we found that the sequences E^^'^^ converge from above and the sequences 
jt;[d,i] fj-om below. We assume that the former provides upper bounds and the latter 
lower ones as in an earlier treatment of the quartic anharmonic oscillator p^[T6] . Thus, 
from sequences of roots with D < 45 we conjecture that 

- 0.698526217161 > ^0 > -0.698526217173, ^^o = ^ 
-2.4713450252412632 > Eq > -2.4713450252412641, vq = 4 
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-0.426405989 > > -0.426405991, t;o = 4 (27) 

We are not aware of any calculation of the eigenvalues and eigenfunctions for this rational 
potential. 

5. Modified RPM for critical parameters 

According to the results of Sec. |3] the appropriate boundary condition at threshold is 
given by equation (fT7|) . Therefore, since the RPM moves a pole of f{x) (a zero of ipix)) 
to infinity it seems reasonable to apply the RPM to the function 

g{x) = -— (28) 

X ip [X) 

We thus have 

i^/(x) - fix)gix) + -gix) = Q{E, x) (29) 

X X 

and 

oo 

g{x)=x^gjX^^ (30) 
i=o 

The first coefficients are 



Qo 


2Qi 
3Qo 




a -^^ 
45 ~ 


4Q2 2Qj llQi 
5Qo 9QI 45 


(31) 


go = Qo 






9i =Qi + 






9'2 = Q2 + 


Y^QoQi + Y^Qo 


(32) 



for s = 1. We apply the RPM exactly in the same way and construct the Hankel 
determinants with the coefficients gj: H'^{E,vo) = \gi^j+ci-i{E , vo)\fj^^ = 0. 

We obtain convergent sequences of roots of H'^{E = 0,vo) = for all the models 
discussed above. In particular. Table [T] shows the first critical parameters for the 
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Gaussian well estimated from the roots of the Hankel determinants with D < 45, d = 
and d = 1. For comparison we add the results of Liverts and Barnea [T3J for the s- 
states of the central-field model. The critical parameters for the central-field model 
with angular momentum quantum number / = are exactly those for the odd states of 
the one-dimensional case. 

6. Central-field models 

The results of Sec. [5] suggest that present approach may also be suitable for the s-states 
of other central-field models. Although the present paper is focused on one-dimensional 
parity-invariant models we can outline a strategy for the treatment of central-field 
models. We write the radial part of the dimensionless Schrodinger equation as 



Qir)- 



ip{r) = 



2r2 

Q(r) = 2[E- V{r)] , ^(0) = (33) 
and assume that 

oo 

Qir) = E Qj^' (34) 
i=-i 

As in earlier papers we define [15] 



r ip[r) 

but in order to apply the RPM to the calculation of critical parameters we resort to the 
auxiliary function 

9{r) = -- ^ 36 
They are related by 

^^?(r) - fir)gir) + lf{r) - Q(r) = (37) 
and can be expanded in a Taylor series about the origin as 



ir 



j=0 



9{r) = T.9jr^ (38) 

j=0 
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As an illustrative example we choose the Yukawa potential 

V{r) = (39) 

and show the results in Table |2] for the first s-states estimated from roots of the Hankel 
determinants with D < 45, d = and d = 1. Present results agree with those of Liverts 
and Barnea |14j and Singh and Varshni [22] up to the last digit reported by them. 

Although present results are more accurate than those of Liverts and Barnea [H] 
and Singh and Varshni [22j one should not conclude that the RPM is superior to the 
approaches developed by those authors. Those other methods are more general and 
apply to all the states of the quantum wells. However, the RPM is a straightforward 
simple approach that applies to a wide variety of problems. In many cases it yields quite 
accurate results and may be suitable for testing other approaches and even for setting 
benchmark data. 

7. Perturbation theory about the threshold 

We can expand the exact energy ( IT9|) for the modified Poschl- Teller potential in a Taylor 
series about f o,n and obtain the perturbation series about the threshold 

{2n + lf {2n + lf {2n + if {2n + if ^ 
^ =Vq- Vo^n (40) 

that converges for all |^| < {2n + 1)^/8. Note that the perturbation correction of first 
order is zero for all states and that we obtain a negative energy for both vq > fo,™ and 
Vq < fo,„ if ^ is sufficiently small, in spite of the fact that n-th state moves into the 
continuum in the latter case. We can carry out a similar calculation for models that are 
not exactly solvable by means of the RPM. In what follows we illustrate the strategy by 
means of the apparently most favourable case of the rational potential f l23|) for which 
the RPM yields the exact solution at threshold. 

The roots of a Hankel determinant Hf){E, vq) = give us approximations to either 
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E{vq) or vq{E). If we substitute 

E = ^(i^e + E^^k^ + ... + E^^\^ (41) 

and vq = vo,n + ^ into the Hankel determinant, then we can obtain the coefficients E^^\ 
j = 1, 2, . . . /c of the perturbation series, the accuracy increasing with D. Based on the 
Hellmann-Feynman theorem 

'^--(^^) (42) 
we expect that 

^(1) = lim ^ < (43) 

for a physically acceptable solution. Since the solutions are not square integrable when 
^ = then the expectation value in equation ( l42l) is meaningless at threshold; however 
the limit fH3|) may hopefully be finite. In fact, E^^^ = for the Poschl- Teller potential. 

It follows from the discussion in the subsection 14.31 that Ei = when vq = 3/2. 
However, if we substitute vq = 3/2 + C, and the series (HTl) into the Hankel determinants 
for s = 1 we obtain the unphysical result 

^ 1 29 ^3 1847 ^4 275357 ^. ^.^e, 

E = H H + 0(^^) (44) 

8 64^ 768 184320 77414400 ^ ^ 

According to this expansion the energy increases as vq increases beyond fo,i = 3/2 in 
contradiction with fH21) and (1431) . This result reflects the fact mentioned above that the 
RPM does not yield the energy Ei for vq close to threshold. 

If we repeat the calculation for the even states we obtain a perturbation expansion 
with the expected slope at threshold: 

8 192^ 11520 19353600 8128512000 ^ ^ 

At first sight, this result is surprising because no new even state should appear when 
3/2 < Vq < A (the ground state remains bound for all fo > 0). The explanation is that 
the RPM yields a spurious even-state energy associated to the divergent solution ipi^Q. 
For example, when vq = 1.51 the RPM with D = 10 and d = yields the actual ground- 
state energy Eq ^ -0.70483 and the spurious root W ^ -0.00124114797000675832 . 
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The considerably greater accuracy of the latter is due to the fact that the RPM yields 
the exact result when vo = 3/2. The question remains why the RPM does not yield the 
energy Ei of the odd state in view of the fact that the calculation of the critical parameter 
is also exact in this case. The perturbation expansion (jH]) with the wrong slope at 
threshold also describes a spurious root. For example, when vq = 1.49 the roots of the 
Hankel determinants with D = 10 and d = yields W ^ -0.0012609753609799139195 
[s = 1) and Eq ^ —0.692231 (s = 0). It is clear that the RPM favours the unphysical 
solutions; in fact, it is also interesting that the spurious roots of the Hankel determinants 
follow the unphysical expansions ( l45l) and f H4l) for both ^ < and ^ > and in either case 
the rate of convergence of the corresponding sequences is remarkably large. For example, 
the expansion fHSj) predicts a positive root for a negative value of ^ and the RPM already 
yields it quite accurately. When vq = 1.49 we obtain W = 0.0012588560223263235359 
in agreement with that perturbation series. 

The second excited state vanishes when vq = 4. The RPM perturbation expansions 
obtained by substitution of f o = 4 + ^ are also unphysical. For example when s = we 
obtain a series with the wrong slope at threshold: 

32^ 2304^ 331776 59719680 1203948748800 ^ 

This results is not surprising if we take into account that the RPM fails to give us the 

second excited state when vq is slightly larger than 4. When s = 1 we obtain the exact 

power series for another unphysical root of the Hankel determinants 

JLf2 1921 3 1186027 4 2551967839 5 o(3^)(A7) 
32^^768^ ^552960^ 696729600^ "^14046068736000^ ^ ^ 

From the roots of the Hankel determinants for vq = 4.01 we obtain Ei ^ —0.429395 
(the actual energy for the first-excited state) and the spurious eigenvalue W ~ 
—0.00031158508464057747545 which is associated to the divergent function 'ip2,i when 

e^o (E^o-). 

For some unknown reason the RPM yields the unphysical roots associated to the 
divergent states more accurately than the physical ones stemming from the convergent 
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states. However, the approach is still a useful tool for obtaining the eigenvalues and 
critical parameters of one-dimensional wells as already shown above. The modified RPM 
discussed in Sec. [5] also yields the same spurious roots; therefore, we may conclude that 
such an unexpected behaviour is inherent in the Hankel determinants constructed from 
either the coefficients fj or Qj. 

Although the RPM yields the whole set of critical values with either the even or 
odd functions for the modified Poschl- Teller, the Hankel determinants do not exhibit 
spurious roots in this case. 

8. Conclusions 

Despite of being a local approximation the RPM may be a useful tool for the calculation 
of critical parameters of one-dimensional and central-field quantum wells. In some 
cases, like the modified Poschl- Teller and rational potentials, the original version of 
the approach yields accurate results. In other cases, like the Gaussian potential, it 
is necessary to resort to a modified algorithm that applies to one-dimensional models 
as well as to the s-states of central-field ones. The accuracy of the estimated critical 
parameters is satisfactory for the Gaussian and Yukawa potentials. 

In the case of the modified Poschl- Teller and rational potentials the RPM yields the 
whole set of critical parameters for both the even and odd solutions to the Schrodinger 
equation. Half of them are associated to the divergent solutions. The occurrence of 
unphysical results is due to the fact that the RPM does not take into account the 
asymptotic behaviour of the eigenf unctions explicitly. For some unknown reason the 
RPM seems to favour the divergent solutions in the case of the rational potential. This 
undesirable behaviour is not a serious limitation because there is no doubt about which 
roots are spurious. 

The Poschl- Teller and rational potentials are different in the sense that the RPM 
yields accurate energies close to the threshold in the former case but not in the latter 
one. Therefore, the three one- dimensional potentials discussed in Sec. H] reveal three 
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different behaviour patterns in the apphcation of the RPM to simple one- dimensional 
parity-invariant quantum wells. 

The RPM yields the exact solutions for the rational potential (!23|) when E = 0. 
This particularly fortunate situation enables one to try perturbation theory about the 
threshold. Surprisingly, the RPM yields exact perturbation series only for the unphysical 
case of divergent functions. The reason may be that the physically meaningful solutions 
do not exhibit power-series expansions about threshold. However, we know that such 
expansions exist in the case of the Poschl- Teller potential as shown in Eq. (HOj) . In this 
case the Hankel determinants do not exhibit spurious roots. 

We have also carried out calculations for potentials of the form V{x) = —vo/{l + 
a;2)™, where m = 5/2,3,4. In these cases the RPM fails to provide the critical 
parameters and the modified RPM exhibits convergent roots. The even and odd critical 
parameters appear separately as in the case of the Gaussian potential discussed in 
subsection 14.21 

9. Appendix 

If we change the independent and dependent variables in the Schrodinger equation with 
the rational potential fl23|) according to 



X 



iz 



(f{z) = ip{iz) = Vl — z'^w{z) 



then we obtain 



d_ 

dz 



d 



This equation is a particular case of the spheroidal differential equation 



d 



w = 



(48) 



(49) 



(50) 



with A = 0, 7^ = —2E and /i^ = 2t>o -|- 1. We can also relate equation f l49l) with the 
associated Legendre equation 



d 



d 
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when E = and u = 0. 

However, for the present discussion we prefer to proceed in a different way. The 
rational potential (123|) exhibits singularities at x = ±i. If we substitute 

^{x) = {l + x"^)" u{x) (52) 

into the Schrodinger equation we obtain an eigenvalue equation for the new dependent 
variable u{x): 

^ 1 + x^ 

2Ex^ + 2 {2a^ - a + e)] m =0 (53) 



+ 



We remove the third term by choosing a to be any one of the roots of 

2a^ - 2a - Wo = (54) 
Then, we expand m in a Taylor series about the origin 

oo 

u{x) = J2 CjX^^^" (55) 

j=0 

and derive a three-term recurrence relation for the coefficients 

(2j + s + l)(2j + s + 2)9+1 + [(2j + s + 2a) (2j + s + 2a - 1) + 2E] 9 
+2Ecj_i = (56) 



1) because of the 



The radius of convergence of this series is unity (limj^oo 
singularities at x = ±i. When E = the recurrence relation becomes a two-term one 
and we can obtain polynomial solutions for particular values of a. Note that if 

{aUk, s) = —k — ^ 
(57) 
a2{k, s) = -k - ^ + 1 

then Cj = for all j > k.We thus have two sets of critical potential parameters: 

(1)/, N _ (2fc+s-l)2-l 



(2a -1)2-1 vyik,s 



\vi'\k,s) = ^2!^ 

(12) 

where Vq = fo(ai,2)- It is interesting to note that the critical parameters exhibit a 
kind of degeneracy: 

vl^\k,l) =vl^\k,0) = ^^^,k = l,2 
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vi,'\k + 1,0)= vj^\k, 1) = + k = l,2,... (59) 

similar to the one found by Joseph P^[2U] for the central-field version of this model. In 
the present case one of the degenerate solutions is convergent (~ 1) and the other one 
is divergent (~ x). The connection between both models becomes apparent if we take 
into account that the states of the central-field model with angular-momentum quantum 
numbers / = — 1 and / = become the even and odd states of the one-dimensional one. 
In order to understand the results derived in the subsection 14.31 by means of the 

RPM simply note that 

^ s _ jj'jx) ^ s _ 2ax _ 
X -ipi^x) X I + x^ u(x) 

is a rational function of x. 
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Critical parameters for quantum wells 

Table 1. Critical parameters for the Gaussian well 

n vo,n (present) (5/2 (Liverts and Barnea) 

1 1.342002325462045769136522184 1.3420023 

2 4.32454875170105636793006 

3 8.89784977356695359409760 8.89785 

4 15.05314025436583553156 

5 22.7867399600521324219 22.78674 

6 32.09666656038554310 

7 42.98170019005867753 42.9817 

8 55.4410239055636496 

9 69.474057353841774 69.47406 

10 85.0803698581927 

11 102.259630867595735 102.25963 

12 121.011579785235 

13 141.336006623056 141.33601 

14 163.23273906943 

15 186.7016335411 186.70163 

16 211.742568809 

17 238.3554413509 238.35544 

18 266.540161872 

19 296.29665268 296.29665 

20 327.62484562 

21 360.52468057 360.52468 



[22] Singh D and Varshni Y P 1984 Phys. Rev. A 29 2895. 
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Table 2. Critical parameters for the Yukawa potential 
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State vo^n (present) (3/2 (Liverts and Barnea) 1/S (Singh and Varshni) 

Is 0.83990388669822801527775556 0.83990390 0.839903886698226 

2s 3.2236301610682666483973 3.2236302 3.22363017 

3s 7.1710139208439285832 7.1710140 7.17101392 

4s 12.6858299220239072 12.685830 12.685830 

5s 19.76942118485633 19.769421 19.769421 

6s 28.42243219866088 28.42243 28.422432 

7s 38.645227430528 38.645227 38.645227 



